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Abstract 

During the past year, a minimal renormalizable supersymmetric SO(10) model has been 
proposed with the following properties: it predicts a naturally stable dark matter and neutrino 
mixing angles # atm and #13 while at the same time accommodating CKM CP violation among 
quarks with no SUSY CP problem. Suppression of proton decay for all allowed values of tan (3 
strongly restricts the flavor structure of the model making it predictive for other processes as 
well. We discuss the following predictions of the model in this paper, e.g. down-type quark 
masses, and neutrino oscillation parameters, U e3 , <5mnsp> which will be tested by long baseline 
experiments such as T2K and subsequent experiments using the neutrino beam from JPARC. 
We also calculate lepton flavor violation and the lepton asymmetry of the Universe in this 
model. 



1 Introduction 



The existence of non-zero neutrino masses appears to have considerably narrowed the choice of 
grand unified theories, with those based on the group SO(10) being preferred for several reasons. 
SO (10) is the minimal group with all the ingredients for a small neutrino mass since its 16- 
dimensional spinor representation contains the right-handed neutrino, vr, which is needed to 
implement the seesaw mechanism^ along with the other fermions of the standard model. It 
also contains the B — L symmetry needed to keep the right-handed neutrino masses below the 
Planck scale and provides a group theoretic rational for the belief that neutrinos most likely 
are Majorana particles. The unification of all the fermions into the 16-plet also raises the hope 
that the number of parameters required to describe fermion masses and mixing will be less than 
those in the standard model, thus making the model predictive for neutrinos. 

Initial attempts to realize this hope were made in a class of minimal SO(10) models with a 
single 10 and a single 126 Higgs multiplets[3El which led to predictions for neutrino mixings as 
well as two of the masses without any extra symmetry assumptions. The gratifying result was 
the natural manner in which the large solar and atmospheric mixings arose in these models with 
predictions in gross agreement with current neutrino observations. The detailed predictions for 
the case where the Yukawa couplings are CP conserving are however away from the current 
central values of the neutrino parameters though still in agreement with observations at three 
u level. 

Encouraged by this initial success in understanding large neutrino mixings, attempts were 
made to study CP violation in this modeljU |3] by making the Yukawa couplings complex. It 
was found that fitting fermion masses and mixings forces the CKM phase to be in the second 
quadrant rather than in the first. This question has been reanalyzed in two recent papers: (1) 
using type II seesaw in these context of SO(10) models^, if one allowed a very high value 
of the strange quark mass, the CKM phase could be in the first quadrant 0; (2) in type I 
seesaw (and mixed case), it is shown that all the masses and mixings (within 99% CL) can 
be fit for a finely tuned range of parameters jH]. While one could take this to be an indication 
possibly of new physics contributions to CP violation, a more conservative point of view would 
be to demand that the model be extended to generate CKM CP violation and see whether its 
predictivity in the neutrino sector is still preserved. One way to achieve this is to go beyond 
this minimal Higgs structure by including the 120 Higgs field [§J EJ. In particular, it was 
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shown in Ref. [HI EI] that the model not only accommodates CKM CP violation but it also 
solves the SUSY CP problem^] and has the potential to solve the strong CP problem as well. 
This would make such models plausible candidates for a theory not only of neutrino masses 
but also of CP violating phenomena. 

As is well known, GUT models generically lead to an unstable proton and therefore the 
lifetime of proton provides additional constraints on them. Several GUT models have already 
been severely constrained by the present experimental bounds on the proton lifetime [T2*| ITS] . 
The question of proton decay in this class of SO (10) models has been under scrutiny in several 
papers[Tn | ll4 | ITH]. In particular, in Ref-^Oj, we showed that requiring the suppression of proton 
decay for both small and large tan j3 severely constrains the textures of the Yukawa couplings. In 
particular, 120 becomes a necessity for this to happen. The reason for this is that cancellation 
of different terms in the decay amplitudes for different proton decay processes are absent in this 
formalism. This makes the model very predictive. In this paper, we study these predictions. 

In the process of suppressing the decay amplitude, the number of free parameters becomes 
less and the model becomes very predictive in the quark-lepton sector. For example, the 
down-type quark masses, \U e3 \, #mnsp etc, are predictions of this model. We also predict the 
lepton asymmetry which then get converted to baryon asymmetry via sphaleron process in the 
thermal leptogenesis scenario. Since the baryon asymmetry is determined quite accurately from 
the recent experimental data, it is important to calculate the prediction of this model. We can 
also predict the probability of muon type neutrino to electron type neutrino oscillation {P u ^ Ve ) 
as a function of distance and the energy of the neutrino beam. The future measurements of 
this probability at the T2K [TH] and the subsequent experiments would shed more lights on 
this model. 

This paper is organized as follows: The model and the Yukawa matrices in our model 
are discussed in section 2, and the natural realization of proton decay suppression and the 
preferable Yukawa structure are discussed in section 3. In section 4, the predictions of our 
model are presented and the impact of these solutions for T2K and lepton flavor violations are 
discussed. The lepton asymmetry is calculated and then it is converted into baryon asymmetry 
in section 5. Section 6 is for the conclusion. 
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2 Fermion Mass Matrices in a Minimal SO (10) Model 

We first introduce the Yukawa interactions and the contents of Higgs fields in the SO (10) 
model. The Yukawa superpotential involves the couplings of 16-dimensional matter spinors ipi 
(i denotes a generation index) with 10 (if), 126 (A), and 120 (D) dimensional Higgs fields: 

1 1 — 1 

W Y = -hij^jH + -fijipitpjA + -h^^D. (1) 

The Yukawa couplings, h and / are symmetric matrices and h! is an anti-symmetric matrix 
due to the SO(10) symmetry. The Higgs doublet fields not only exist in H, A, D, but also 
exist in other Higgs fields which are needed in the model. For example, a 210 Higgs field (<&) 
is employed to break the SO(10) symmetry down to the standard model and $ contains Higgs 
doublet fields. One 126 Higgs multiplet A is introduced as a vector-like pair of A and this field 
also contains a Higgs doublet. The VEV of this pair reduces the rank of SO(10) group and helps 
to keep supersymmetry unbroken down to the weak scale. Altogether, we have six pairs of Higgs 
doublets: (p d = (Hj°, D l dl D d ,~A d , A d , <p u = (if* , D*, D%, A u , A u , $ u ), where superscripts 
1, 2 of D Ut d stand for SU(4) singlet and adjoint pieces under the G422 = SU(4)xSU(2)xSU(2) 
decomposition. The mass term of the Higgs doublets is given as (<Pd)a(MD)ab{<Pu)b, and the 
expression of the matrix Md is given in R.ef. [17] . The mass matrix of the Higgs doublets is 
diagonalized by unitary matrices U and V: UMdV t = M^, iag . We assume that the lightest 
Higgs pair (MSSM doublets) has masses of the order of the weak scale. The MSSM Higgs 
doublets are given as linear combinations: Ha = Ui a ((pd) a , H u = V* a (ip u ) a . Since we concentrate 
on the structure of Yukawa couplings, we do not describe the dynamical reason of the mass 
hierarchy in this paper. 

We use "Y-diagonal basis" (or SU(5) basis) to describe the standard model decomposition 
of the SO(10) representation[r7| ITS] . The above Yukawa interaction includes mass terms of the 
quark and lepton fields as follows: 

jymass = hHf(qd c + le c ) + hHl°(qu c + lv c ) (2) 
+ -^fMqd c - 3ie c ) + ^/A„(gn c - 3tv c ) + V2fv c v c A R + V2f££A L 

+ UD\{qd c + ie c ) + h'Dl(qu c + iu c ) + -^ti D 2 d (qd c - 3£e c ) - -=ti ' D 2 u (qu c - 3£u c ), 

v 3 v 3 

where q, u c , d c , £, e c , v c are the quark and lepton fields for the standard model, which are all 

unified into one spinor representation of SO(10). The VEVs of the fields Ar : (1,1,0) and 
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Al '■ (1,3,1) give neutrino Majorana masses. We obtain the Yukawa coupling matrices for 
fermions as 

Y u = h + r 2 f + r 3 h', (3) 
Y d = nCh + f + h'), (4) 
Ye = r l (h-3f + c e h'), (5) 
Y u = h-3r 2 f + c v h', (6) 

where the subscripts u, d, e, v denotes for up-type quark, down-type quark, charged-lepton, and 
Dirac neutrino Yukawa couplings, respectively, and 

h = V n h, f = U 14 /(V3n)f, h! = (U 12 + U 13 /V3)/ ri ti, (7) 
U u V l5 V 12 -V 13 /V3 

n = ^ r2 = ri ^' r3 = r V 12 + f W3' (8) 

U 12 - V3U 13 V 12 + V3V 13 fn . 

Ce = 7=i C v = T\ = . (9) 

+ U 13 /V3 u 12 + u 13 /Vs 

The light neutrino mass is obtained as 

m|j ght = Ml — M?Mr 1 (M?) t , (10) 

where M? = Y U (H U ), M L = 2V2f(A L ), and M R = 2^2f(A R ). 

In this paper, we consider that the Yukawa coupling matrices for quarks and leptons are 
hermitian. The hermiticity of fermion mass matrices can be obtained in the assumption: The 
original Lagrangian has symmetry under charge conjugation in addition to the SO(10) gauge 
symmetry. In the Y-diagonal basis, the SO(10) vector representation in terms of X a can 
be decomposed such that (Xx, X 3 , X 5 , X 7 , X s ) transforms as 5-plet, and (X 2 , X 4 , X 6 , X 8 , X ) 
transforms as 5-plet under SU(5)xU(l) decomposition. We define the charge conjugation of 
the field written in the Y-diagonal basis, X a «-> X*. The conjugation of the higher rank tensor 
representation is also defined similarly using the F-diagonal basis. Since the SO(10) symmetric 
Lagrangian has parity invariance as an internal symmetry (actually, it is D-parity), the theory 
now has CP symmetry. Imposing that Lagrangian is invariant under the CP conjugation, the 
Yukawa couplings, hy, and hL and all masses and couplings in the Higgs superpotential 
are all real. However, when SO(10) symmetry is broken down, the CP symmetry can be 
spontaneously broken by the VEV of 45 Higgs field. Due to the non-existence of cubic terms for 
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45 Higgs field, the VEV of 45 can be pure imaginary. Consequently, the mixing of the lightest 
Higgs doublets with the Higgs doublets present in 120 involves a pure imaginary coefficient 
which will make the fermion masses hermitian in this modeljU]. This symmetry has wider 
implications. For instance, the model presents a solution to the SUSY CP problem and strong 
CP problem. It is possible to explain the quark masses and mixing angles and the neutrino 
sectors by using the above parameters 

3 Proton Decay and Flavor Structure 

As mentioned above, it was shown recently by us^Oj that the suppression of proton decay for 
all tan/? determines the flavor structure of the three matrices h, f, h' to a very narrow range. 
We review this argument in this section. 

The proton decay is mediated by the colored Higgs triplets, <pt + '■ 1> — 1/3) + c.c.) 
and (p c + <Pc '■ (( 3 > !» ~ 4 / 3 ) + c - c -)- Tnese Hi gS s triplets appear in 10+120+126+126+210 
multiplets. We generate both LLLL (Cl) and RRRR (Cr) operators: 

-W B = \cf l q kqiq ^ + Cf l etu^. (11) 

These operators are obtained by integrating out the triplet Higgs fields, ip^ = (-£%, -Dy, D^, Ay, 
Ay, A^, $y) and y? T = (H T , D T , D' T , Ay, A T , A T , $ r ). The fields with '" are decuplet, and the 
others are sextet or 15-plet under SU(4) decomposition. The Cr operator has also contributions 
from the triplets, ip^ = (D^, A^r) and ip c = (-De? Ac). The mass term of the Higgs triplets are 
given as ((Pj^) a (M T ) ab (ip T ) b + (<Pq) a (Mc) abi^c) &■ The mass matrices, My and M c , are 7x7 and 
2x2 matrices respectively, and their explicit forms are given in the literature [IZ]. The Yukawa 
couplings which cause proton decay are written as 

Wy ip - = hH T (q£ + u c d c ) + hH T (^qq + e c u c ) + /Ay (q£ - u c d c ) + f~E T (^qq - e c u c ) 
+ V2fA' T e c u c + V2h'{D T u c d c + D' T qi- D T e c u c + D' T e c u c ) 
+ 2fA c d c e c + 2h'Dc u c u c + 2h'D c d c e c . ( 12) 

The dimension five operators are written by the Yukawa couplings h, f and h! as follow: 

C % l kl = chijhu + Xifijfki + x 2 hijfti + x 3 fijh k i + xji'^h^i + x^h'^f^i, (13) 
C% kl = c hijh k i + yifijfki + y2hijfki + yzfijh k i + y^Kj^ki + ysKjfki 

+ yehijh'u + yrfijKi + VzKjKi + VsKifjk + yioKih' jk . (14) 

5 



The coefficient c is given as c = (M^ 1 )u, and the other coefficients Xj, y^ are also given by the 
components of M^ 1 or Mq 1 . Note that the and Ay have opposite D-parity and we get 
y 3 = —x 3 . The y 9 and y w terms are generated by <pc + fc- 

The proton decay operators can be written conveniently by diagonalizing the Higgs triplet 
mass matrix Mt by two unitary matrices, X and Y, as XMtY t = diag(Mi, M 2 , ■ • • , M7), 

C L kl = E TF^ h + X ^f + V2X a3 h%(Y al h + Y a5 f) kh (15) 

°t l = E Tr( X *i h ~ X -^' + V^X^/iOy (^ai/i - (^a5 " v 7 ^)/ + V2(Y a3 - Y a2 )h') kl 

+ (2/9, 2/10 terms). (16) 

In the SU(5) limit, only one of the colored triplets is much lighter than the others, i.e. Mi -C M a 
(a 7^ 1), and we can obtain the following relations for the diagonalizing matrices from the 
explicit form of the Higgs mass matrices in Ref . [T71 ITH] : Un = X n , Vu = Y u , C/ 14 = X 14 = 0, 
V 15 :Y 15 :Y W = V3:1:V2, U 12 : U 13 : X 12 : X 13 = V 12 : \\, : Y 12 : Y 13 1 : : v '2 : v '2. 
Using the above relations, we get 

r 2 -> 00, / - 0, r 3 = 0, c e = -1, c = 2nFi 2 /C/ 12 , (17) 

for the Yukawa matrices in Eqs.(j3J|BJ) and thus, as expected, we get the SU(5) relations, Y u = Yj, 
Yd = Y e T , and the dimension five proton decay operators can be written in terms of the Yukawa 
couplings as C % [ kl ~ C^ M ~ (Y d )ij(Y u ) k i/Mi and we know that proton decay is not suppressed 
at this limit and this limit is also not good to satisfy the quark-lepton masses. 

The proton decay amplitude can be written as A = a 2 f3 p /(4:7rMTmsusY)A, where 

A = cA hh + XiAff + x 2 A hf + x 3 A fh + x 4 A h , h + x 5 A h , f . (18) 

The coefficients c and Xi are given in Eq. (jl3|) . and there are also similar Cr contributions. To 
satisfy the current nucleon decay bounds, we need ^ 10~ 8 , I^U-^I <, 2 ■ 10~ 8 and 

|^4 n ^icp| ^ 5 ■ 10 -8 if the colored Higgsino mass is 2 • 10 16 GeV, and squark and wino masses are 
around 1 TeV and 250 GeV, respectively. Instead of inducing any cancellation among different 
terms in the amplitude, we will try to suppress the individual contributions. 

One way to suppress the decay amplitude is by demanding cancellation among different 
terms, a strategy common in the literature. In order to achieve that, we need a cancellation 
among h, f and h' to have small couplings for first and second generations in the expressions in 
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Eqs. (J13ll4|) . However, we also need cancellation among the same couplings to generate the large 
mass hierarchy among the quark masses and in general, the coefficients r 2 , r 3 in up- type Yukawa 
matrix and in proton decay operators are unrelated. Further, the 126 Higgs contribution 
has opposite signatures (7/3 = — £3) for Cl and Cr. Therefore, the cancellation required to 
obtain small Yukawa coupling for Y u by tuning r 2 / can not simultaneously suppress both Cl 
and Cr operators by tuning Xu in Eqs. (jl51llfi|) . Moreover, the 120 contribution to the proton 
decay amplitude has vanishing contribution to the kl part of C % [ kl due to antisymmetric h! . 
Thus, if the cancellation in Y u requires a tuning of r 3 h', the Cl operator can not be suppressed. 

The best way to avoid the cancellation is to choose smaller values of r 2j 3. Let us start with 
r 2i3 ~ 0. In this case, the Cl can be written as C l [ kl oc (Y u + r yh f ) i j(Y u )ki and in the operator 
C l ^ kl , kl part is also related to Y u . This will correspond to the case where X14, Y\$ ~ 0. The 
RRRR contribution to p — > Ku T mode is suppressed compared to the minimal SU(5) model by 
a suppression factor X u /Xd ~ 1/100 for tan/5 ~ 50. Similarly, since the kl part of Cl are also 
related to the Y u instead of Yd, the LLLL contribution to the p — > Kv is also suppressed even 
for tan/3 ~ 50, compared to the SU(5) model (since A c /A s ~ 1/5). However, these suppressions 
are not enough to satisfy the current experimental bound. 

In order to satisfy the bounds naturally, we need A hh <, 5 • 10~ 8 in the expression, Eq.(|18|). 
If <^ 10~ 7 , we need to tune x^ and t/i for every decay mode to cancel Ahh, which is 
unnatural. The A^h depends on the magnitudes of the elements from the [1,2] block of h which 
is specifically determined from the fit to the up-type Yukawa coupling as a function of r 2 and r 3 . 
In the case r 2 .3 ~ 0(1), the suppression of up- and charm-quark Yukawa couplings are acquired 
by fine-tuning of r 2) 3, and thus the elements from [1,2] block of h are of the order of down- and 
strange-quark Yukawa couplings. In that case, we find Ahh ~ 10~ 4 which requires a very high 
level of fine-tuning for all the decay modes. In the case r 2 , r% ~ for generical fits, the Yukawa 
coupling h is close to Y u in the ^-diagonal basis by definition. However, even in the case, 
the Ahh is of the order of 10~ 7 for generical fits and further tuning among the coefficients xi 
and yi is needed to satisfy the current experimental data. We therefore need a specific type of 
Yukawa texture to suppress the proton decay rate. To suppress Ahh, the elements hn and /i 22 
(in /i-diagonal basis) are needed to be suppressed rather than the up- and charm-quark Yukawa 
couplings, respectively. As a result, we need Yukawa texture to be h ~ diag(~ 0, ~ 0,O(l)). 
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Once h is fixed, the other matrices / and h' are almost determined as 

/ ~ ~ A 3 \ / 



f 








V A 3 



~ A 3 \ 
A 2 A 2 
A 2 A 2 / 



h' 



-A 3 
-A 3 



A 3 A 3 \ 
A 2 
-A 2 ) 



(19) 



where A ~ 0.2. The correct charm mass is generated by r2m s /mb ~A C ( [r*2] —0.1 — 0.15), and 
down-quark mass and Cabibbo angle 9c are generated by h' 12 with m^/m s ~ sin 2 9c- Then we 
need /n & 0(A 6 ), /12 ^ 0(A 4 ) and r 3 ^ 0(A 2 ) to obtain proper size of up-quark and electron 
masses. We have r\ ~ m^/rrit tan /?. In the basis where F u is diagonal, Ahh in this texture is 
not completely zero but can become much smaller than 10~ 8 . 

We show one example for numerical fit for tan/3(M^) = 50 : h = diag(0, 0, 0.638), 



f 



( -0.00044 0.00208 \ 

-0.00044 0.00945 0.0101 
V 0.00208 0.0101 0.0071 / 



h! 



f -0.0022 0.00046 \ 

0.0022 0.0181 

V -0.00046 -0.0181 / 



(20) 



n = 0.966, r 2 = 0.135, r 3 = 0, |c e | = 0.987. 

After we suppress the Ahh, we also need to examine the contribution of the other components 
e.g. Aff t hfjh,h'f,h'h,---- Their coefficients, Xi,yi, involve the colored Higgs mixings, which can be 
suppressed by our choice of the vacuum expectation values and the Higgs couplings. According 
to our numerical studies, some of the mixing angles must be less than about a few percent in 
the case of tan (3 ~ 50 to suppress the decay. However, the mixing angles can become larger as 
tan/3 becomes smaller. In the above example of numerical fit, p — >• Kv^ and p — > Kv T modes 
are dominant for LLLL operator, and for RRRR operator, respectively. The Ahh for p — > Kv^ 
mode is ~ 2 • 10~ n . The amplitudes for other components are 10~ 8 — 10~ 6 . 

We have seen that the suppression of proton decay requires r 3 ~ 0, which is same as 
SU(5) condition for 120 Higgs coupling, and r 2 — 0.1 — 0.15, which however is not a SU(5) 
condition. The second condition, as well as suppression of colored Higgs mixing, is implemented 
by requiring Uu ^> V15, A 14 , Y 15 . We note that if U ~ V is satisfied (which is like a up-down 
symmetry and therefore tan/3 ~ 50), the r 3 ~ condition also satisfies the SU(5)-like condition 
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(21) 



for 120 Higgs coupling in Y e and Y u , though these conditions are not required to be satisfied 
to fit the fermion masses or to suppress the proton decay. 
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Figure 1: \U e s\ is plotted as a function of m s . The two regions are described in the text. 

4 The Model Predictions 

Now we are ready to discuss the predictions of the model. The number of parameters in the 
models is 17: 3(h), 6(f), 3(h') and 5 Higgs parameters (7*1,2,3, c e>u ). (We are working in the 
diagonal h basis). Now in the process of explaining the proton decay, some of the parameters 
are redundant to fit masses and mixings. We choose ^11,22 = and r% = 0. The r% can be 
0(X 2 ) in our Yukawa texture, but such a value of r 3 gives only a small correction to the CKM 
mixings. We use a very small fn, but even in the case fn = 0, the obtained Yukawa texture 
can fit all masses and mixings. We make it a free parameter since it affects the leptonic fitting. 
Since we will be working type II seesaw (actually, the Yukawa texture of / is compatible to 
bi-maximal mixings), c„ is redundant in fitting fermion masses and mixings. This reduces the 
number of parameters to 13. In order to fix the remaining 13 parameters, we use the up-type 
quark masses, charged lepton masses, the CKM angles and the phase, the ratio of the squared 
of neutrino mass differences (Am^/Am^), and the bi-maximal mixings as input parameters. 
Consequently, the down-type quark masses, U e3 and <5mnsp etc are the predictions of this model. 
We itemize below the essential features of our predictions. 
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Figure 2: \U e s\ is plotted as a function of V u b and Am^ ol / Am\. 
4.1 Strange Quark Mass 

In order to see how the down quark masses and in particular the strange quark mass gets 
constrained in our model, we first note from Eqs. (|4l5ll9J) we get det M e ~ |c e | 2 det ., when 
fn <C A 5 and thus, the down quark and the charged lepton mass matrix leads to the Georgi- 
Jarskog (GJ) relation naturally when the SU(5) relation (c e = —1) is satisfied in 120 Higgs 
mixing. 

The strange quark mass is predicted since no cancellation happens to derive V c b (due to 
r 2,^3 *C 1 and the hermiticity condition). The predicted value of strange quark mass has two 
separate regions, roughly m s ~ 1/3 m^(l ± 0(A 2 )). In our numerical calculations, the negative 
sign corresponds to a strange quark mass: m s (/x = 2 GeV) ~ 120 — 130 MeV. (The strange quark 
mass at 1 GeV is obtained by multiplying the strange mass at 2 GeV with 1.35.) This strange 
mass value is in the range of lattice derived value, m s (/x = 2 GeV) = (105 ± 25) MeV[IH]. The 
GJ relation is realized for this result and thus the fitting of the lattice result supports |c e | ~ 1 
(when fn <C A 5 ). If we use the positive signature in the rough estimation, we find the following 
value of the strange quark mass, m s (/i = 2 GeV) ~ 155 — 165 MeV which is allowed by the 
QCD calculation^] ( in this case, non-zero values of fn are needed to make |c e | = 1). In 
Figure n we plot \U e $\ against the strange quark mass. We see the two regions corresponding 
to two different values of strange quark masses (corresponding to two different signs). We also 
note that the larger values of strange mass prefers lower values of U e s- We use a s (Mz) = 0.118 
and a 3-loop QCD and 1-loop QED running to evolve the strange mass from 2 GeV scale to 
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the electroweak scale) 

Our prediction of strange-down quark mass ratio is m s /m^ ~ 17 — 18, 19 — 20.5. The 
non-lattice value of the ratio is m^/m^ = 18.9 ± 0.8|21j. 

The m u /md ratio is not a prediction of this model since m u is an input. To obtain a large 
atmospheric mixing, ^33 needs to be suppressed, and thus bottom-tau Yukawa coupling needs 
to be unified within several percent. In order to achieve such a situation, large tan/5 ~ 50 is 
needed (or tan/3 ~ 2). Thus, the bottom mass is predicted. However, this prediction is lost 
since it gets a sizable contribution from the soft SUSY breaking terms for large tan f3. 

4.2 \U e2 \ 

Since there is no cancellation, we get the following stable approximate relation for f/ e3 : 

We,? - , ta ?J?L Am s 2 ol / Am 2 A . (22) 



We also have the following relation since U e % is related to the ratio /13 and /23 : 

\TJ I ~ 1 Vub 



V 



cb 



(23) 



These relations are obtained since /12 is small compared to /13. Actually, when / 12 — /13, U e -s 
can be canceled to be zero without relating to the solar mixing angle and V u b- In our proton 
decay suppressed texture, the fu is suppressed due to the assumption hn <C X u , and thus the 
\U e s\ prediction can be stable. 

The \U e s\ is bounded from below through the mass squared ratio and solar mixing angle 
and has an upper limit due to the experimental bound on \V u b\ = (3.67 ± 0.47) x 10~ 3 . In 
FigureEl we plot \U e s\ as a function Vub and Am^ ol / Am\. We find that for \V u b\ < 0.0045 and 
AmlJAm\ > 0.022, the \U e3 \ is bounded to be 0.06 - 0.11. We plot l^sl as a function \c e \ 
in Figure El The value of c e = —1 (corresponding to the SU(5) condition for 120) gives rise to 
\U e s\ ~ 0.1. Both large and small strange quark masses are allowed for c e = —1. 

4.3 MNSP Phase 

The MNSP phase is also stable in this model and is given by the approximate expression: 

1 sin^f, / -1 c e h\ \ , , 

smW ~7rf si T a " szrj' (24) 
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Figure 3: \U e ^\ is plotted as a function of c e . 

where 9\ 2 and 9^ are mixing angles in the diagonalizing matrix of Y e and neutrino mass matrix, 
respectively. Approximately, sin #^3 ~ \U e z\. We plot \U e ^\ as a function of MNSP phase <5mnsp 
in Figure |3J We find that the 5mnsp lies i n the 2nd or 4th quadrant if c e < 0. The 1st and 
the 3rd quadrants are absent since the solar mixing angle becomes small due to a cancellation 
between Q\ 2 and Q\ 2 . In Figured we plot |c e | as a function of <5mnsp- We see that c e = — 1 
(SU(5) relation of 120 Higgs mixing, and also with a favorable value by the GJ relation), 
generates sui^mnsp — ±(0.5 — 0.7). To obtain this result, we assume that type II contribution 
is dominant. 

The location of <5mnsp in the 2nd or 4th quadrant has impact on the probability of to v e 
oscillation {P v -^ Ve ) which will be measured at the T2K experiment and at the newly proposed 
Tokai-to-Korea experiment |22j. This probability depends on sine and cosine of <5mnsp> distance 
(L), energy of the neutrino beam, mass squared differences (Am 2 ^, AmJ 2 ), 3 mixing angles, 
and matter density. We take the energy of the beam is 0.7 GeV and the values of mass squared 
differences: Am 2 3 = 2.5 x 10~ 3 eV 2 , Am 2 2 = 8 x 10~ 5 eV 2 and U e3 = 0.1. In Figure we 
plot the probability as a function of distance. The probability for S = 330° is about 1.8 times 
bigger compared to the probability for 5 = 135° when the beam arrives at Kamioka from Tokai 
(L = 295 km). The difference is magnified much more if we have a detector installed at Korea 
(L = 1000 km). Also we notice that a peak in the distribution appears at T2K and at the 
Tokai-to-Korea experiment. The location of the peak however will change if we change the 
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Figure 4: \U e3 \ and c e is shown as a function of <5mnsp f° r c e < 0- 



above parameters. 

5 Leptogenesis 

Let us consider the calculation of baryon asymmetry of the Universe. As we will see be- 
low, a lepton asymmetry is generated by CP violating decays of right-handed neutrinos out 
of equilibrium |23j. The sphaleron processes |24j convert the lepton asymmetry to a baryon 
asymmetry |25|. Our goal is to see if the recent experimental value of the baryon-to-photon 
ratio, i.e. [2^] 

7] B = (6.3 ±0.3) x 10~ 10 , (25) 

can be reproduced in our model. 

The baryon-to-photon ratio is obtained as 

Vb = -aj^Kd, (26) 

where a is a sphaleron conversion factor a ~ 0.35, and Iff® is coming from the fact that the 
photon number density increases as the degrees of freedom present at the epoch of leptogenesis 
annihilate and is given by /° = g*/g°, where g* is the number of degrees of freedom at T = 
(mass of the lightest right-handed neutrino) and g° = 3.91, giving 1//° — 0.016. The e± is the 
amount of CP violating lepton asymmetry via the lightest right-handed neutrino, which is 
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L (km) 



Figure 5: The probability P v ^ Ue function of L(km). 

defined as 

e v r(iv 1 -> mi) - r(iv 1 -» 4gj 

The k is an efficiency factor, which does not related on the CP violation of lepton asymmetry 
and parameterizes the effects of scattering and decay processes. In the thermal leptogenesis 
scenario, the k is a calculable number by solving Boltzman equation. The k is a function of an 
effective neutrino mass, fhi, which defined as 

*, = MMi, (28) 

M Nl K J 

where is the Dirac neutrino mass matrix in the basis where right-handed Majorana mass 
is diagonal with real and positive eigenvalues, Mjv r The Ref. |27j shows that the k for 10~ 2 eV 
< rhi < 10 3 eV is approximately given as 

.^^Vh^r 8 , (29) 



mi J V 10 3 eV 
and in Ref.j2H], the k is given by a power law as 



, ,'0.01 eV\ L1±01 . , 

k = (2±1)x10" 2 — = . (30) 



mi 

Now let us see the prediction of our model. As we have seen, flavor structure has been 
already fixed in our model, and thus the ei and k can be calculated as a function of the 
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lightest right-handed neutrino mass and SU(2)^ triplet Higgs mass. For example, since the 
\(Mj?y M®]n is almost fixed in our model, the effective neutrino mass is a function of the 
lightest right-handed neutrino, M Nl . As a result, the efficiency factor k depends on only M Nl . 
Consequently, the recent experimental data, Eq. ()25j) . gives us a prediction to the mass scale. 

Let us discuss the possible contributions to the CP violating lepton asymmetry in our model. 
The decay amplitude of the lightest right-handed neutrino in one-loop involves the right-handed 
neutrinos and SU(2)x triplet Higgs [28] present in 126 + 126 and 54 (if there is one). The 
right-handed neutrino loop contribution is obtained as 

-± ^^^A F (m, (31) 



where Y v is the Dirac neutrino Yukawa coupling in the basis where right-handed Majorana mass 
is diagonal with real and positive eigenvalues, M^ i . Since the expression of the does not 
depend on the overall scale of right-handed Majorana neutrino mass, the value of is almost 
determined in our model. The value is proportional to c u , which is a coefficient of 120 Higgs 
contribution to Dirac neutrino Yukawa coupling, Eq.©. Our prediction is [2E] 

|ef| = (l-4)xl0" 4 (c„/2). (32) 

Suppose that this dominate the lepton asymmetry ei, the value of k is determined to be k ~ 
(0.25 — 1) x 10~ 3 (2/c u ) to satisfy the current experimental data. The values of k correspond to 
the effective neutrino mass, rhi = (0.1 — 0.4) {c u /2) eV, and we obtain the lightest right-handed 
Majorana mass is M Nl = (0.4 — 1) x 10 13 (2/c v ) GeV. Using the relations, / = V3ri/Uuf, 
M R = 2v / 2/-y R , we obtain the corresponding value of VEV is v R ~ Uu (2/c^) (1 — 2.5) x 10 15 
GeV. 

We also have the triplet Higgs loop contribution, ef , and the contribution may disturb the 
above prediction of lightest right-handed neutrino mass. In fact, as discussed in Ref. |28t l3*U]. 
this triplet Higgs loop contribution may dominate t\ naively when we consider type II seesaw. 
However, the amount of lepton asymmetry depends on the flavor structure, and such naive 
estimation may not be correct. Actually, in our model, the can be calculated as a function 
of lightest right-handed Majorana mass, and the triplet mass. To see this, let us see how we 
can realize type II seesaw. 

In order to satisfy the type II VEV magnitude for v L = (A L ), one 54 Higgs is needed |3*T|. 
There are then two pairs of SU(5) submultiplet 15s in the theory (15 Higgs of SU(5) is the 
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one that contains the triplet that leads to the type II seesaw). For the triplet VEV term to 
dominate, one of these two pairs must have a mass around a scale of 10 13 GeV or so. It was 
shown in Ref.[3*T] that a linear combination of the SU(5) 15 sub-multiplets in 54 and 126 can 
have a mass around 10 13 GeV without conflicting with the symmetry breaking and coupling 
unification and the other one becomes closer to the SO(10) breaking scale. 

We can write down the superpotential involving the triplets and their interactions (for 
simplicity, we omit the terms from 210 Higgs): 

W = \ X H»H U & E + A 2 ^A S A L + X 2 v R A E A L + Y A LLA L + M A A L A L + M E A E A E , (33) 

where Y A = 2\/2/, and A^ is a triplet in 54 Higgs. The VEV is given as 

v 2 

v L ~\ lSA c A —^, (34) 
where s A , c A are mixing (sin and cos) of SU(2)l triplets, 

M A 2 

We denote M Al2 as eigenmasses of SU(2) L triplets and we assume M Al -C M A2 to make type 
II dominant and M A to be around GUT scale. Note that M Al is a free parameter even if we 
fix the v E . However, assuming that the Higgs couplings are less than 0(1), the M Al should be 
less than around 10 13 GeV. 

The lepton asymmetry via triplets in the loop is given in the Ref. |281 l3*Uj . 



3 x M Nl lm(YjY A Y t 



ef = --^A lCAgA 7; vi A1 g{ v ) (36) 
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3 M Nl lm(Yjm l jY;) n 



Svr vl (YjY u ) n 

where G(y) = y In and y = M\J M^^. Our prediction of the triplet loop contribution is 

14 1 = (2-5 - 6, x 10- M 2) (^) (^) G (V) . (37) 

We note that the function G(y) is always smaller than 1. If the triplet mass M Al is 10% 
compared to the lightest right-handed Majorana mass, Mjv i: we obtain G{y) ~ 0.046 and 
the triplet loop contribution can be negligible compared to the right-handed neutrino loop 
contribution e^. When M Al and are comparable, then the both triplet and right-handed 
neutrino loops contribute to the lepton asymmetry. 



16 



Q_ 



to 
c 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 




Figure 6: sin 5 M nsp is plotted as a function of the ratio of the type I contribution to the type 
II contribution (|(m*) 33 1/(^)33) in the heaviest neutrino mass [Left]. The |c e | is plotted as a 
function of 5mnsp for r] B > and [Right] for v R ~ U 1A (2/c v ) (1 - 2.5) x 10 15 GeV. 



It is possible to check the predictions, Eqs. (|32I37|) . using the Yukawa couplings shown as an 
example in section 3. the Dirac neutrino coupling in the above basis with c v = 2 (which is a 
favorable value as we noted in the last paragraph in section 3) is given by: 

/ 0.002 0.003 exp(-1.54i) 0.0026 exp(-0.344z) \ 
Y u = -0.0167 0.021 exp(-1.53z) 0.025 exp(3.37z) • (38) 
V -0.229 0.417 exp(4.70i) 0.422 exp(-3.019i) / 

The Majorana neutrino coupling matrix which gives the right-handed neutrino masses is pro- 
portional to the / coupling, 



/ = diag. (0.00108, 0.00303, 0.0185) . 



(39) 



If there is no huge cancellation between and ef\ the resulting lepton asymmetry E\ needs 
to be O(10 -4 ) and consequently, the efficiency factor k is O(10~ 3 ) and the lightest right-handed 
Majorana mass, M^, is 10 13 GeV. Note that we also have solutions where the Majorana mass 
scale is larger. In this case, since k becomes larger, the lepton asymmetry e\ needs to be 
canceled between and by choosing the triplet mass. 

We note that the triplet Higgs decay also can produce lepton asymmetry by its own decay, 
but this contribution vanishes when the Hermitian structure of Yukawa coupling is assumed. 

When Mtv x ~ 10 13 GeV, the type I contribution can disturb the prediction of oscillation 
parameters obtained in the previous section. We find that the \U e ?\ prediction is not changed 
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much, but the prediction of the MNSP phase is modified due to the phase of type I contribution. 
As shown in figure |H1 the sin <5mnsp splits from the type II prediction when type I contribution 
increases (However, the correct neutrino fit in our model can not be obtained by pure type 
I contribution). In figure El we show the sukSmnsp as a function of the ratio of the type I 
contribution to the type II contribution (1(^)331/(^)33) to the heaviest neutrino mass. We 
use the h, f and h! values presented in section 3 (we use c e ~ —1 and c v = 2) and vary vr/U^ 
continuously. In figure El we also show the MNSP phase for the case of positive r\R dominated 
by the contribution given in Eq.([31j). When vr/Uu > 2 x 10 16 GeV, type I contribution can 
be negligible and the MNSP phase lies around 135° — 150°, 315° — 330° (when c e ~ —1) as we 
have seen in the previous section. In this scenario of large vr, the generated lepton asymmetry 
may be an order of magnitude bigger compared to the observation since the efficiency factor 
k becomes larger. However, the left-handed triplet Ax, contribution in the loop can produce a 
10% level cancellation and generate the right amount of lepton asymmetry. The location of the 
MNSP phase, which can be observed at T2K and the subsequent experiments as mentioned 
before, can be a probe to distinguish the scenario. 

In our discussion so far, we did not concern ourselves with the nature of SUSY breaking 
mechanism. As is well known, gravitino production puts an upper bound of about 10 9 GeV on 
the reheating temperature Tr after inflation when the gravitino mass is in the range, 100 GeV 
^ m 3/2 ^ 1 TeV[32]. Since in our model, the lightest right-handed neutrino which is supposed 
to be responsible for leptogenesis has a mass of 10 13 GeV or so, clearly the thermal leptogenesis 
picture outlined in this section will not work in this case. The possibilities to make the thermal 
leptogenesis scenario available are to consider a light gravitino such as 7713/2 < 16 eV[33], late 
time entropy production in a gauge mediation model[33], quasi thermal picturel^S]. Otherwise, 
one can consider a non-thermal leptogenesis where the inflaton plays a key rolej^H]- We do not 
give details of the implications of such a choice for our SO (10) model except to note that it 
does not affect the fermion sector of the model. 

6 Lepton Flavor Violation 

We now discuss the lepton flavor violating processes e.g. \i — > ej, r — > fi'j etc. The operator 
for li — > lj + 7 is: 



£^. 7 = — lj a^ v q v (aiP L + a r P R ) k ■ + h.c. (40) 
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Figure 7: BR[/i — > e + 7] and BR[r — > \i + 7] are plotted as a function of mi/2 for different 
values of m . 

where -Pl ; ,r = (1 =F 7s)/2 and cr M1/ = | [7^, 7"]. The decay width for Zj — > Zj + 7 can be written 
as: 

r(z^/ J+ 7) = ^(hr + kl 2 ). (4i) 

Then the branching ratio is obtained by multiplying this decay width with the lifetime of the 
U lepton. The supersymmetric contributions include the neutralino and chargino diagrams. 

We work in the basis where the charged lepton masses are diagonal at the highest scale of 
the theory. We use the mSUGRA boundary conditions, i.e. m : universal scalar mass, m 1( / 2 : 
universal gaugino mass, A : universal trilinear mass term, tan j3, sign of fi to calculate the 
BRs. Since we have At, ~ A T , the masses of the Higgs are not unified with the universal scalar 
mass for the sfermions for realistic parameter space after including the finite SUSY one loop 
correction to the b-quark mass. We plot BR[/i — > e + 7] and BR[r — » fi + 7] as a function of 
mi/2 for different values of m in Figure and we find that the BR[/i — > e + 7] can be quite 
large. 
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7 Conclusion 



The suppression of proton decay is a major challenge for the grand unified models. We have 
constructed a minimal SO(10) model which can suppress the proton decay naturally. The model 
has 10, 126 and 120 Higgs multiplets to generate the fermion masses. The CP symmetry of 
the model keeps the fermion masses hermitian at the grand unified scale and the SUSY CP 
problem is also under control. The neutrino masses are generated by type II seesaw. The model 
has 13 parameters to fit quark and lepton masses and mixings, and thus it gives predictions 
for down-type quark masses, \U e s\ and 5mnsp- The prediction for |{7 e3 | is 0.06 — 0.11 with the 
upper limit being imposed by V u b and the lower limit by the mass squared ratio (Am s 2 ol /Am^) 
and the solar mixing angle. The strange quark mass can be small (lattice calculated size) 
or large (QCD calculated size). The phase 5mnsp is i n the 2nd or 4th quadrant when the 
SU(5)-like condition is satisfied in the 120 coupling. This feature has distinctive impact on the 
probability of to v e oscillation [P v ^ Ve ) which will be measured at the T2K and subsequent 
experiments. The lepton asymmetry generated in the decay of the right-handed neutrinos in 
this model produces the observed amount of baryon asymmetry. Depending on the magnitude 
of the symmetry breaking scale Vr, the observed baryon asymmetry predict the 5mnsp phase 
in different quadrants. The lepton flavor violating Br \p — > e + 7] can be large in this model. 
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